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Abstract 

This paper addresses a structural design problem in control systems, and explicitly takes into 
consideration the possible application to large-scale systems. More precisely, we aim to determine and 
characterize the minimum number of manipulated state variables ensuring structural controllability of 
switched linear continuous-time systems. Towards this goal, we provide a new necessary and sufficient 
condition that leverages both graph-theoretic and algebraic properties required to ensure feasibility of 
the solutions. With this new condition, we show that a solution can be determined by an efficient 
procedure, i.e., polynomial in the number of state variables. In addition, we also discuss the switching 
signal properties that ensure structural controllability and the computational complexity of determining 
these sequences. In particular, we show that determining the minimum number of modes that a switching 
signal requires to ensure structural controllability is NP-hard. 


I. Introduction 

Switched systems have been intensively studied, and the primary motivation comes partly 
from the fact that these systems have numerous applications in the control of mechanical 
systems, process control, automotive industry, power systems, aircraft and traffic control, and 
many other fields [13], [27]. Among switched systems, those with all subsystems described 
by linear differential equations, referred to as switched linear systems, have attracted most of 
the attention [13]. Recent efforts aimed to analyze controllability and reachability properties 
of these systems [2], [11], [27], [28]. Nonetheless, only recently controllability was studied 
for the class of uncertain switched linear system, i.e., the parameters of subsystems’ state 
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matrices are either unknown or zero [14]. This assumption eopes with seenarios where the 
system parameters are diffieult to identify and obtained with a eertain approximation error. 
Thus, struetural properties that are independent of a speeifie value of unknown parameters 
are of partieular interest. Subsequently, a switehed linear system is said to be strueturally 
eontrollable if one ean find a set of values for the unknown parameters sueh that the eorresponding 
switehed linear system is eontrollable in the elassieal sense [14]. Whereas in [14] neeessary and 
suffieient eonditions to eharaeterize the struetural eontrollability were provided, due to eeonomie 
eonstraints we aim to determine and eharaeterize the smallest subset of aetuated state variables 
yielding struetural eontrollability. 

The problem of determining the smallest subset of state variables yielding eontrollability is 
eommonly referred to as minimum controllability problem, and it was studied in the eontext 
of linear time-invariant systems in [15], [25]. Reeently, the minimum eontrollability problem 
has also been explored under additional energy eonstraints, in partieular, metries that depend 
on the eontrollability Grammian, see [3], [4], [12], [17], [26], [29]. Alternatively, if we aim to 
aehieve struetural eontrollability, we refer to the problem as structural minimum eontrollability 
problem. This problem has been fully addressed in the eontext of linear-time invariant systems 
in [19] for homogeneous eosts; and the eomputational eomplexity analyzed for several elasses 
of systems in [1]. In the eurrent manuseript, we extend these results to the ease of switehed 
linear eontinuous-time systems. Notiee that whereas in [1], [19] graph-theoretie properties of 
struetural eontrollability rely in direeted graphs interpretations of the system, these no longer hold 
to eharaeterize struetural eontrollability of struetural switehed linear eontinuous-time systems, 
see [14] for details. In partieular, the analysis of struetural eontrollability of struetural switehed 
linear eontinuous-time systems eannot be redueed to the analysis of a struetured linear system. 
Consequently, in this paper we provide a systematie approaeh that leverages the eombination of 
graph-theoretie and algebraie eonditions to obtain and eharaeterize the solutions to the struetural 
minimal eontrollability problem for struetural switehed linear eontinuous-time system. 

The eurrent work also differs from [23], [24], where the struetural minimum eontrollability 
problem aimed to ensure struetural eontrollability for eaeh mode of the switehed eontinuous- 
time linear system; note that this eonservative notion eontrasts with the eontrollability definition 
eonsidered in the present manuseript. In [20] the struetural minimal eontrollability problem 
for linear time-invariant systems was eonsidered under heterogenous eost, i.e., the variables 
aetuated ean ineur in different eosts. In partieular, in [20] this problem is shown to be poly- 
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nomially solvable, and in [16], the eomputational eomplexity was improved when binary eosts 
are eonsidered. More reeently, in [22] the problem was extended to the case where a state 
variable has a cost that depends on the input that actuates it, hence, leading to a multiple 
heterogenous cost scenario. Alternatively, the problem of determining the minimum number of 
actuators from a given collection of possible actuator-state configurations was shown to be (in 
general) NP-hard [18]. Notwithstanding, in [21] it was shown that the same problem can be 
polynomially solvable when the dynamic matrix is irreducible. 

The main contributions of this paper are fourfold: (i) we provide a new necessary and sufficient 
condition that leverages both graph-theoretic and algebraic properties required to ensure structural 
controllability of switching linear continuous-time systems; (ii) we characterize the solutions to 
the structural minimum controllability problem for switched linear continuous-time systems. In 
particular, we characterize dedicated solutions, i.e., an actuator can only actuate a single state 
variable, and minimal solutions, i.e., the minimum number of actuators actuating the minimum 
number of state variables; (Hi) we propose an algorithm that leverages both graph-theoretic 
and algebraic properties of structural controllability of switching linear continuous-time systems 
to determine a solution in (mn)“, where n denotes the number of state variables, m denotes 
the number of modes of the switching linear continuous-time system, and a < 2.373 is the 
exponent of the n x n matrix multiplication; and (iv) we provide new insights on how the 
switching sequences affect the controllability of a structural switching linear continuous-time 
system. In particular, we show that determining the minimum collection of nodes in a sequence 
of modes ensuring structural controllability is NP-hard. 

The rest of the paper is organized as follows: Section II provides the formal statement of 
the problem addressed in this paper. Next, Section III reviews some concepts, introduces key 
results in structural systems theory and establishes their relations to graph-theoretic constructs. 
In Section IV, we present the main results. Next, we present an illustrative example in Section V. 
Finally, Section VI concludes the paper and discusses avenues for further research. 

II. Problem Statement 

In this section, we formally introduce the structural minimum controllability problem for 
switched linear continuous-time systems. 
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Consider the following switehed linear eontinuous-time system 

x{t) = (1) 

where a : IR+ —)■ M = {1,. .. ,m} is a switehing signal, x{t) G MA the state of the system 
at the instanee of time t, and u{t) G MA represents the pieeewise eontinuous input signal. In 
the sequel, we identify (1) by the pair that eontains m modes with subsystems 

{Ai, Bi), i G {1,..., m}, and a{t) = i implies that the ith. subsystem (Aj, Bi) is aetive at time 
instanee t. Further, the switehed linear eontinuous-time system (1) is said to be eontrollable 
(or equivalently, {A^(t), Ba-(t)) is eontrollable) if for any initial state a;(0) = xq, and a desired 
state Xd, there exists a time instanee tf > 0, a switehing signal a : [0,f/) —)■ M and an input 
u : [0,f/) —)■ sueh that x{tf) = Xd- This notion of eontrollability enables the analysis of 
switehing systems where we either have aeeess to ‘eommon’ transitions and knowledge of the 
existing modes of the switehing system, or the eases where the eontroller is equipped with 
supervisory eapabilities enabling the system to switeh between modes. 

Due to eeonomie eonstraints, one is interested in deploying the minimum aetuation eapabilities 
that enable the eontrollability of the system, whieh ean be eaptured by the following optimization 
problem. Given the switehed linear eontinuous-time system (1), we aim to determine the sparsest 
input matriees {B*}'^^ for the different m modes required to ensure eontrollability, as a solution 
to the following optimization problem: 

m 

min XI ll^illo 

s.t. (A^(fc),5^(fc)) is eontrollable, 

where ||M||o is the zero (quasi) norm, i.e., it eounts the number of non-zeros entries in matrix 
M. Unfortunately, this problem is NP-hard even when m = 1, see [15] for details. 

Furthermore, the parameters assoeiated with the linear time-invariant modes are often aeeu- 
rately known, whieh motivates the use of struetural system theory [8]. Struetural linear systems 
are linear parameterized systems with a given strueture, i.e., the entries of the state spaee matrix 
are either free parameters or fixed zeros. Let A^-^t) G {0,1}"^^” denote the zero/nonzero strueture 
or structural pattern of the system matrix whereas B„(t) ^ {0,is the struetural 

pattern of the input matrix B^p)- More preeisely, an entry in these matriees is zero if the 
eorresponding entry in the system matriees is equal to zero, and deseribed by an arbitrary 
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parameter (denoted by one) otherwise. Therefore, a pair B^(^t)) is said to be structurally 

controllable if there exists a pair respecting the structure of B^(^t)), i-C-, 

same locations of zeros and nonzeros, such that is controllable. Further, it can be 

shown that if a pair {A^(^t), B^(^t)) is structurally controllable, then almost all (with respect to 
the Lebesgue measure) pairs with the same structure as {Aait)-, B„(t)) are controllable [14]. In 
essence, structural controllability is a property of the structure of the pair {Aa-(t), Bcr(t)) and not 
of the specific numerical values. 

Subsequently, the structural minimum controllability problem for switched linear continuous¬ 
time systems problem can be stated as follows: 


Vi Given the structure of the matrices of the switched linear system in (1), i.e., {Ai}'i[Li, we aim 
to determine the sparsest collection of input matrices required to ensure its structural 

controllability, i.e., that is the solution to the following problem: 


min 


s.t. 


E mio 


2 = 1 

(yio-(fc), i?o-(fc)) is struct, controllable. 


Notice that a solution to Vi may consist of Bi with columns with all entries are equal to zero, 
which can be disregarded when considering the deployment of the inputs required to actuate the 
system. In addition, in the worst case scenario, we obtain structural controllability by taking the 
identity matrix as the input matrix, which justifies the dimensions chosen for the solution search 
space. Furthermore, some solutions may comprise at most one nonzero entry in each column; 
in other words, solutions in which each input actuates at most one state variable. These inputs 
are referred to as dedicated inputs, and they correspond to the columns of the input matrix Bi 
with exactly one nonzero entry. Additionally, if a solution is such that all its nonzero 

columns consist of exactly one nonzero entry, it is referred to as a dedicated solution, otherwise 
it is referred to as a non-dedicated solution. 

Finally, note that the solution procedure for Vi also addresses the corresponding structural 
observability output matrix design problem by invoking the duality between observability and 
controllability in linear time-invariant systems for each mode of the switching linear continuous¬ 
time system [10]. 
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III. PRELIMINARIES AND TERMINOLOGY 

In this section, we review some notions of controllability of switched linear continuous-time 
systems, and their counterpart using structural systems theory [8]. 

To assess the controllability for switched linear continuous-time systems consider the following 
definitions. 

Definition 1 ([14]). The controllability matrix for switched linear continuous-time system as 
described in (1) is given by 

C{Aa-{k)j = [-Bi, B2, ..., Bm, AiBi, A2B1, ..., 

AraBl, A1B2, A2B2, . . . , AmB 2 , . . . , AiBm, 

^2-Bm; • • • ; A^Bm, AfBi, A2A1B1, . . . , A^AiBi, 

A1A2B1, A2B1, . . . , • • • , AiAmBm, 

A 2 A^B^, ..., AlB^,A^^B,, A 2 Ar^B^, 

..., ..., A^^B^, 

4, 4^-2 TJ 4”“1 R 1 

. . . , J->mi ■ ■ ■ 1 ^rn 


O 


Additionally, we have the following result. 

Theorem 1 ([14]). The system described by (1) is controllable if and only if rank C(Ao-(fc), B^^k)) = 
n. o 

Now, we associate with the pair (Aj, Bi), with Aj, Bi G (0, a directed graph (digraph) 

V{Ai, Bi) = {Vi, Si), referred to as the system digraph, with vertex set Vi and edge set Si, where 
Vi = UiU Xi with Xi = {x\,... ,xU and Ui = {u\,... ,uf} represents the state and input 
vertices, respectively. In addition. Si = Sx,,Xi A Su,,x, where Sxi.x, = {{xi,xi) : [Ai]jk A 0} 
and Sui,Xi = • [Bi]jk AO} represents the state edges and input edges, respectively. 

Similarly, we can define a state digraph V{Ai) = {Xi,Sxi,Xi)- A directed path is a sequence of 
directed edges where every edge ends in a vertex that is the starting of another edge and no vertex 
is used twice. A state vertex is said to be non-accessible by an input vertex if there exists no 
directed path from an input to the state vertex. Later, given matrices Mi,..., Mm G {0, 


June 10, 2016 


DRAFT 



7 


their structure can be combined into a single matrix M in terms of M = Mi V ... V Mm, 
where V corresponds to the entry-wise operation where if at least one of the entries is non-zero, 
then it provides a non-zero entry, and zero otherwise. In addition, [Mi,..., M^] denotes the 
concatenation of matrices Mi,..., Mm- 

Next, we introduce the notion of a bipartite graph associated with a mi x m 2 matrix M given 
by B{M) = {C,TZ,£cji), where IZ = {ri,... ,rmi} and C = {ci,... , 0 ^ 2 } correspond to the 
labeling row vertices and column vertices, respectively, and Sc,^ = : Mij ^ 0}. The 

bipartite graph is an undirected graph with vertex set given by the union of the partition sets C 
and TZ, which we refer to as left and right vertex sets, respectively. A matching M. C £c,n is 
collection of edges that have no vertex in common. A maximum matching is a matching with 
maximum cardinality among all possible matchings. For ease of reference, if a vertex in the left 
and right vertex set does not belong to an edge in a maximum matching, we then refer to it a right- 
and left-unmatched vertex, respectively. Additionally, we can consider weights associated with 
the edges in a bipartite graph to obtain a weighted bipartite graph {B{M) = {C,TZ,£c,'R.),w), 
where w : Ec^n —)■ M. Subsequently, we can consider the problem of determining the maximum 
matching with the minimum sum of the weights, that we refer to as the minimum weight 
maximum matching. The minimum weight maximum matching can be generally solvable in 
Cl(max{mi, 7772 }“), where a < 2.373 is the exponent of the n x n matrix multiplication [9]. 

In addition, a digraph Vs = (V 5 ,£^ 5 ) is a subgraph of V = (V,£^) if V 5 C V and 8s E 
8. Finally, a strongly connected component (SCC) is a maximal subgraph (there is no other 
subgraph, containing it, with the same property) Vs = (V 5 , £s) of V such that for every u,v eVs 
there exists a path from u to v and from v to u. We can create a directed acyclic graph (DAG) 
by visualizing each SCC as a virtual node, where there is a directed edge between vertices 
belonging to two SCCs if and only if there exists a directed edge connecting the corresponding 
SCCs in the digraph V = {V,8), the original digraph. The DAG associated with V{A) can be 
computed efficiently in 0{\V\ + |£^|) [6]. The SCCs in the DAG may be further categorized as 
follows: an SCC is non-top linked if it has no incoming edge to its vertices from the vertices 
of another SCC. 

Finally, consider a mi x m 2 matrix M, and let [M] = {P G . p.. = g if Mij = 0}, 

then the generic rank (g-rank) of M is given by g-rank(M) = max rank(P). 

P&[M] 

Now, we revisit necessary and sufficient conditions for the structural controllability of switched 
linear continuous-time systems. 
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Theorem 2 ([14]). A switched linear continuous-time system (1) is structurally controllable if 
and only if the following two conditions hold: 

(i) V{Ai V ... V Am, i?i V ... V Bm) has no non-accessible state vertex; 

(ii) g-rank {[A^,Am, i?i, • • •, Bm]) = n. o 

We note that whereas verifying if the conditions in Theorem 2 hold can be done efficiently [14], 
designing the sequence of sparsest input matrices such that those conditions yield is a more 
challenging problem. In fact, a greedy strategy may consist in sequently try to ensure each 
condition. Nonetheless, optimality of such strategies cannot (in general) be ensured. Therefore, 
one should resort to such strategies only when the problem at hand is computationally intractable, 
for instance, NP-hard. In this paper, we will show that Vi can be polynomially solvable by 
leveraging both graph-theoretic and algebraic characterizations of the conditions in Theorem 2 
captured by the following results. 

Lemma 1 ([7]). Let M G {0, exists a maximum matching of B{M) with size n 

if and only if g-rank{M) = n. o 

Lemma 2 ([19]). The digraph V{Ai V ... V Am, i?i V ... V Bm) has no non-accessible state 
vertex if and only if there exits an edge to a state vertex in each non-top linked SCC of the DAG 
representation of V{Ai V ... V Am) from an input in V{Ai V ... V Am, i?i V ... V Bm)- o 

IV. Main Results 

In this section, we present the main results of this paper. More specifically, we characterize 
all the solutions to Vi. This goal is achieved in three steps. First, we determine a dedicated 
solution that enables structural controllability by performing actuation into a single mode, i.e., 
we determine a dedicated solution B* such that B^ = B and B* = 0 for i = 2,. ■ ■ ,m ensures 
structural controllability of {Aa-(t), B*^^-^) (see Algorithm 1 whose correctness is provided in 
Theorem 4). Second, we describe in Theorem 5 the non-dedicated solutions Bf and B* = 0 
for i = 2,... ,m, which can be obtained from the dedicated solutions B*. In the last step, in 
Theorem 6, we consider the former characterization to describe all possible solutions to Vi- 
Finally, we present a discussion of results regarding the switching signals required to ensure 
structural controllability of {A„(t), B*^^^-^). In particular, we show that determining the minimum 
sequence of modes in such switching is NP-hard. 
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We start by leveraging graph-theoretic and algebraic conditions presented in Lemma 1 and 
Lemma 2, to rewrite Theorem 2 as follows. 

Theorem 3. A switched linear continuous-time system (1) is structurally controllable if and only 
if the following two conditions hold: 

(i) there exits an edge to a state vertex in each non-top linked SCC of the DAG representation 
of V{Ai V ... V Am) from an input in V{Ai V ... V Am-, i?i V ... V Bm); 

(ii) there exists a maximum matching of 

B{[Ai, ..., Am, 5i, ..., Bm\) with size n. o 

The conditions provided in Theorem 3 can now be used to obtain a dedicated solution to 
Vi- More specifically, we propose Algorithm 1 to obtain B such that Bi = B and Bi = 0 for 
all i = 2,... ,m is a dedicated solution to Vi. Towards this goal, Algorithm 1 consists of a 
two-step algorithm that determines the smallest collection of state variables required to ensure 
both conditions in Theorem 3. In the first step, it find the set containing the maximum collection 
of state variables that simultaneously contribute to satisfy both conditions in Theorem 3. We 
show that this set can be obtained by considering a MWMM on a weighted bipartite graph. 
The weighted bipartite graph requires the careful crafting of a ‘surrogate’ matrix S that will 
encode the graph-theoretic properties required to ensure condition (i) in Theorem 3, and a weight 
function enables the connection of condition (i) and the algebraic condition in Theorem 3- 
(ii). Consequently, a MWMM determines the maximum set of state variables that need to be 
actuated, while satisfying both conditions in Theorem 3. More specifically. Algorithm 1 finds 
a MWMM of a bipartite graph i3([Ai,..., Am, >§]), where the matrix S has as many columns 
as the number of non-top linked SCCs in the DAG representation of 'D{Ai V ... V Am)- The 
non-zero entries in column i oi S correspond to the indices of the state variables that belong to 
the Ath non-top linked SCC. In addition, we consider weights in the edges of the bipartite graph: 
the edges associated with the nonzero entries of Ai have zero weight and the edges associated 
with nonzero entries in S have unitary weight. In particular, if an edge in the MWMM contains 
the vertex corresponding to the columns of S, then the row vertex r* in the edge indicates that 
the state variable Xi contributes to satisfy simultaneously both conditions in Theorem 3. Hence, 
in the second step, one remains to determine (independently) the smallest sub collection of state 
variables fulfilling both conditions in Theorem 3. 

The next result establishes the correctness and analyzes the implementation complexity of 
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ALGORITHM 1: _ 

Input: Structural linear time-invariant dynamies in eaeh mode of the struetural switehing 
linear eontinuous-time system deseribed by 

Output: Input matrix deseribing a dedieated solution ©(JT”), where 0(^7”) represents the 
n X n diagonal matrix with entries in J' different from zero 

Step 1. Determine the non-top linked SCCs i G X = {1, • • • , (3}, of 
V{Ai V ... V Am) = (A", Sx,x)- 

Step 2. Consider a weighted bipartite graph i5([Ai,..., Am, >S]) = (C, IZ, Scji)^ where S 
is a n X matrix and S'jj = 1 if Xj G AfJ, and the eolumn vertiees be re-labeled as 
follows: the eolumns of Ai are indexed by {c\,... and the eolumns of S are indexed 

by {si,..., S/j}. In addition, let the weight of the edges e G IJ {c\,..., c^} x 7^ be 

equal to zero, and the weight on the edges e G {si,..., xTlhe equal to one. 

Step 3. Let Ai' be the maximum matehing ineurring in the minimum eost of the weighted 
bipartite graph presented in Step 2. 

Step 4. Take J' = {i : {sj, ri) G Ai',j G {1,..., /3}}, i.e., the row vertiees assoeiated with 
S that belong to the edges in the MWMM Ai' (i.e., those with weight one). In addition, let 
J" = {I,... ,n}\{i e {I,... ,n] -. {c],ri) G Ai',k G (1,.. .,m},j G (1,... ,n}}, and 
J'"' eontains the index of a single state variable from each non-top linked SCC with 
pe{l,...,/3}\J'. 

Step 5. Set J = J'vj J" vj J'". 


Algorithm 1. 

Theorem 4. Algorithm 1 is correct, i.e., it provides a dedicated solution {.B*}™!, with 
obtained using Algorithm 1 and B* = 0 (i = 2,... ,m), is a solution to Vi. Furthermore, its 
computational complexity is 0{{mnA I3)A, where a < 2.373 is the exponent of the nxn matrix 
multiplication. o 

Proof: The eorreetness of Algorithm 1 follows from the faet that the indiees in J' identify 
the minimum set of dedieated inputs, simultaneously maximizing the inerease in the g-rank of 
[Ai,..., Am, D(770] with respeet to [Ai, ..., Am] by |77'|, and the dedieated inputs assigned to 
state variables in different non-top linked SCCs. This follows from observing that (by eonstrue- 
tion) B{[Ai, ..., Am]) results in a minimum weight maximum matehing Ai with zero weight 
and size \Ai\', henee, by Lemma 1, it follows that g-rank([Ai, ..., Am]) = \Ai\. Subsequently, 
a minimum weight maximum matehing Ai' of B{[Ai, ..., Am, *§]) equals \Ai'\ — \ Ai\', henee, 
inereasing by \Ai'\ — \ Ai\ the g-rank of [Ai, ..., Am, D(77')] with respeet to [Ai, ..., Am], and 
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contributing to satisfy both conditions in Theorem 3. Nevertheless, it may be insufficient to ensure 
condition (ii) in Theorem 3, which is fulfilled by taking into account the minimum (additional) 
collection of dedicated inputs indexed by J”. In addition, by construction of S it follows that 
D(j7') corresponds to dedicated inputs that are assigned to state variables in different non-top 
linked SCCs. Hence, \J'\ non-top linked SCCs have incoming edges from different inputs in 
the system digraph. Thus, contributing to satisfy condition (i) in Theorem 3, but may not be 
enough to ensure this condition, which is accounted for by considering the minimum collection 
of dedicated inputs indexed by J”' that ensures condition (i) in Theorem 3. 

In summary, the total number of additional dedicated inputs required such that g-rank([Ai, 

..., Ami U J'')]) = n is minimized by considering Step 4. Similarly, the total number of 
additional dedicated inputs 1’^'” required such that there exist no non-accessible state vertices 
in V{Ai V ... V Arm'0{J' U J'")) is minimized by considering Step 4. Notice that 
are not assigned to non-top linked SCCs previously assigned, otherwise they would have been 
considered in ©(J”'). Finally, by setting J = J' U J” U J'”, as in Step 5, (D(J"), 0,..., 0) is 
by construction a solution to Vi, since both conditions in Theorem 3 hold. Further, it is minimal 
since |j7’| is minimal, which implies that ll-®*llo = \J\- 

The computational complexity follows from noticing that Step 2 can be solved using the 
Hungarian algorithm that finds a MWMM of B{\Ai ,..., Am, *S]) in Cl(max{|C|, |7^|}"), where 
\C\ denotes the number of column vertices in B{[Ai ^..., Am, >S]), and a < 2.373 is the exponent 
of the n X n matrix multiplication, whereas all other steps have linear complexity; hence. Step 2 
dominates the final computational complexity, leading to the final complexity of C>(|C|“), since 
|C| > |77|. ■ 

Remark 1. We notice that if the structural switching linear continuous-time system only possesses 
one mode, then it boils down to a structural linear time-invariant, and the characterizations 
obtained in [19] can be retrieved. Contrarily to the approach presented in [19] that is moti¬ 
vated by the graph-theoretic characterization of the system digraph. Algorithm 1 requires both 
graph-theoretic and algebraic characterizations, since graph-theoretic properties for structural 
switching linear continuous-time system are quite diverse from those known for structural linear 
time-invariant, see [14] for details. o 

Next, we characterize all the possible sparsest matrices that are solutions to Vi when a single 
node is actuated. 
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Theorem 5. Given J” and J”' as in Algorithm 1, then Bl = 0(^7”' U J'”) V 0(^7'"), where 
0(X) is the nxn matrix with exactly one non-zero entry in row indexed in J and zeros otherwise, 
and B* = 0 for i = 2,... ,m is a solution to Vi. o 

Proof: The proof follows by noticing that both conditions of Theorem 3 hold. Condition (ii) 
in Theorem 3 holds because a maximum matching of U J”)]), is also 

a maximum matching of B{[Ai,... ,Am-,E\)- Therefore, since the maximum matching of the 
former has size n, the latter also has size n. Secondly, Theorem 3-(i) also holds for 'D{Ai V 
... V Am, B); more precisely, V{Ai V ... V Am, 0(7"' U J” U T”''')) satisfies Theorem 3-(i), see 
Theorem 4, which implies that there exists a directed path from an input to every state variable. 
Therefore, by considering B = D(7^' U J") V 0(7^"') there exist an input edge from an input 
to the same state variables indexed by J'”. Nevertheless, the inputs from where the input edges 
start are indexed by the columns with non-zero entries in 0(7^"'). ■ 

Note that in Theorem 5, the matrix 0(7^'") has no constraints on the number of non-zero 
entries in each column. Hence, B = D(7^' U J”) V 0(7^"') is not necessarily a dedicated 
solution. Subsequently, we can characterize the minimal solutions to Vi as follows. 

Corollary 1. Given J', J” and J'” as in Algorithm 1, then B\ = D(7^'U7^") V M(7^'"), where 
M(X) is the nxn matrix with exactly one non-zero entry in row indexed in J in some column 
indexed by J' U J” and zeros otherwise and B* = 0 for i = 2,... ,m is a minimal solution to 
Vi. o 

Proof: The proof follows from Theorem 5, and noticing that the entries associated with 
columns indexed by J' U J” cannot be shared by the same column, as consequence of the 
construction in Theorem 4; in particular, it would compromise condition (ii) in Theorem 3. ■ 

Finally, we provide the most general characterization of the sparsest input matrices that are 
solution to Vi. 

Theorem 6. Let J', J" and J"' as in Algorithm 1, and B = D(7^' U J”) V 0(7^"'), where 0(X) 
is the nxn matrix with exactly one non-zero entry in row indexed in J and zeros otherwise. 
If B*, for i = 1,..., m, is such that the following holds: 

• B* contains as columns only the columns of B, but only once the non-zero columns; 

• all non-zero columns of B are present in [B ^,..., Blf; 
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• no two Bj, Bl, with j, k = 1,... ,m and j ^ k, contain the same non-zero column vectors 
ofB, 

then is a solution to Vi. o 

Proof: By noticing that under the mentioned conditions the minimality is ensured, i.e., we have 
YlT=i ll-®i*llo = ll-B||o> we only need to show the feasibility of the solution. Towards this goal, ob¬ 
serve that there exists a 2mnx2mn permutation matrix P such that [Ai,..., Am, Bl,. ■ ■, = 

[Ai,..., Am, B, Onxn, ■ ■ ■, Onxn], whcrc Onxn Is the n X n zero matrix. Hence, by invoking 
Theorem 5 both conditions of Theorem 3 hold, and the result follows. ■ 

The main results also provide new insights about the importance of the sequence of modes in 
the transition required to ensure structural controllability of the switching linear continuous-time 
systems. More specifically, we now do the following observations: 

Remark 2. Given a switching signal that ensures structural controllability of the switching linear 
continuous-time systems, the order of the modes among which the systems transitions does not 
impact the structural controllability of it. In fact, this follows from the conditions presented in 
Theorem 3, since both conditions are invariant to permutation. o 

Subsequently, one may wonder which modes are the crucial to ensure both conditions in 
Theorem 3. This problem can be partially understood from the a solution obtained in Algorithm 1, 
which is captured in the following remark. 

Remark 3. Given the minimum weighted maximum matching M' obtained in Step 3 in Al¬ 
gorithm 1, if an edge (c*,.) G AA' then it follows that mode i is being considered as part of 
the switching signal to ensure structural controllability of the switching linear continuous-time 
systems, since it contributes to ensure condition (ii) in Theorem 3. Nonetheless, this is not the 
same to say that there is no other minimum weighted maximum matching AA" where (c*,.) ^ AA". 
In other words, there might exist different switching signals ensuring structural controllability of 
the switching linear continuous-time systems, and these can be partially captured by the minimum 
weighted maximum matchings. Furthermore, it is possible to characterize which edges belong to 
any maximum matching [5], which implies that those modes need to be part of any sequence of 
modes among which the system transitions. Therefore, these modes should be considered as part 
of the design proposed in Theorem 3. In other words, the designer should consider to deploy 
actuation capabilities in the modes that are strictly required in a sequence to ensure structural 
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controllability, i.e., at these modes, the input matrices should be non-zero. o 

Finally, we notice that determining the minimum sequence of modes a switching signal among 
which the system should transition to ensure structural controllability is NP-hard. Formally, 
consider the following problem. 

V 2 Given a structurally controllable determine the minimum number of modes 

m' that a switching signal a{t) should consider to attain structural controllability. o 

Theorem 7. Problem V 2 is NP-hard. o 

Proof: Consider the well known NP-hard problem, the set covering problem, that can be stated 
as follows: given a universe of elements U = {!,... ,n} and a collection of subsets 
with Sj C U, determine the subcollection that contains U, where X C {!,... ,m} and 

there is no other X' such that \X'\ < |X| satisfying the same conditions. 

In order to show that V 2 is NP-hard, we need to polynomially reduce the set covering problem 
to V 2 (see [6] for an introduction on the topic). As a consequence, a solution to V 2 enables the 
reconstruction of a solution to the set covering problem, which implies that finding a solution 
to V 2 is at least as difficult as finding a solution to the set covering problem. Towards this 
goal, we associate with each mode of the switching system a subset Si, and we assume that 
any mode suffices to ensure condition (ii) in Theorem 3; more specifically, we assume that 
Ai e {0, i].("^+i)x(»^+i) and diag(Aj) = [1,..., l](n-i-i)xi for alH = 1,... ,m. Furthermore, let 
A = Ai V ... V Am be such that T>{A) = (V = {xi,..., Xn+i}) is a directed tree rooted in xp, 
hence, by considering = 1 and = 0 for the remaining i,j = 1,... ,n + 1, and 

= 0 for /c = 2 ,..., m, i.e., only the first mode is actuated and a dedicated input actuates xi, 
it follows that condition (i) in Theorem 3 yields. Notice that there are n directed edges in the 
directed tree, which we can enumerate as {ei}r=i- Therefore, each of the m modes can exhibit 
in its digraph representation the edges {ej}j(zs^ for k = 1,... ,m, besides the self-loops in all 
state variables. 

Under the present construction, it is not difficult to realize that any solution to V 2 consists 
in finding the smallest subcollection of modes such that the system digraph V{A,B), where 
A = \li^x*Ai and X* denotes the indices of the modes that contain a directed spanning tree 
rooted in an input. Therefore, it follows that there exists a collection of edges producing such 
tree, which implies that there exist a collection of sets {iSjjiei* covers lA, and the result 
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follows. 


V. Illustrative Example 

Consider a switched linear continuous-time system with three modes where Ai G 

{0,and [^i]i ,2 = [^ 2 ] 3,2 = [^ 3 ] 4,4 = 1 and zero otherwise. Then, we have 


A = Ai V A 2 V As 


0 10 0 
0 0 0 0 
0 10 0 
0 0 0 1 


Hereafter, we aim to determining the sparsest configuration of inputs that renders the system 

structurally controllable, i.e., a solution to Vi. To this end, we consider Algorithm 1 (whose 

correctness and computational complexity are provided in Theorem 4). 

First, the DAG representation of the state digraph associated with 'D{A) contains four SCCs 
depicted by dashed gray boxes in Figure l-(d); in particular, and A/J are two non-top linked 
see. Subsequently, we have S as follows: 


0 0 
1 0 
0 0 
0 1 



Fig. 1. In this figure, we depict the state digraph associated with Ai, with i = 1, 2, 3, in (a)-(c), respectively. The state digraph 
of the union A = Ai V A2 V A3 is depicted in (d) and contains three SCCs depicted by dashed gray boxes; in particular, 
and M 2 are two non-top linked SCC. In (e)-(f), we represent the bipartite graph obtained in Algorithm 1 and the MWMM. 
Alternatively, in (g)-(h) we illustrate the bipartite graph obtained in Algorithm 1 and the MWMM if a single mode is considered 
with A as its dynamics. 


In addition, the weights associated with the edges in A 2 , S]) are as follows: the 

edges that contain the vertices C 13 and C 14 have unitary weight (depicted by the blue edges in 
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Figure l-(e)), and all other edges ineur in zero weight (depieted by the blaek edges in Figure 1- 
(e)). At Step 3 in Algorithm 1, the MWMM is represented by the eolleetion of red edges in 
Figure l-(f), whieh we denote by M\. Lastly, taking Ml, we obtain J"' = {2}, J'" = 0 (i.e., 
there is no need to inerease the size of the maximum matehing), and J'”' = {4} that ensures 
the all nodes in V{A, B){J' U J" U J'")) are aeeessible. 

Consequently, we obtain that (i?i, B 2 , - 83 ) = {p{J' U J” U J'"), 0, 0) is a dedieated solution 
to Vi, by invoking Theorem 4. From Corollary 1, we obtain that (Bl, B'A = {p{J' U 
J”) U M(j7’"'), 0, 0) is a minimal solution to Vi. If we want the aetuation to be distributed 
among different modes of the switehed linear eontinuous-time system, one just needs to reeall 
Theorem 6 ; in partieular, {B”,B'AB'l) = (D(J"'),D(J"'"),D(J"")) is a solution to Vi. 

Remark 4. The solution obtained for structural linear continuous-time switching systems cannot 
be retrieved from that of a structural linear time-invariant systems when a joint state digraph 
of A is considered. More specifically, consider a switching system with a single mode whose 
dynamics is A, then by executing Algorithm 1 , one needs to consider the bipartite graph in 
Figure l-(g), where the edges in blue have unitary weight and the remaining have zero weight. 
Then, a possible MWMM determined in Step 3 in Algorithm 1 contains the edges depicted in 
red in Figure l-(h), which implies that J' = {2}, J” = {3} and J”' = {4}. In other words, an 
additional state variable is required to be actuated, i.e., X 3 . o 

VI. Conclusions and Further Research 

In this brief paper, we provides a new neeessary and suffieient eondition that leverages 
both graph-theoretie and algebraie properties required to ensure structural controllability of 
switching linear continuous-time systems. With this condition we characterize the solutions to the 
structural minimum controllability problem for switched linear continuous-time systems. Further, 
we provided an efficient algorithm that determines a solution to the problem. Finally, we provides 
new insights on how the switching sequences affect the controllability of a structural switching 
linear continuous-time system. In particular, we show that determining the minimum collection of 
nodes in a sequence of modes ensuring structural controllability is NP-hard. Future research will 
focus on considering different actuation cost per state variables, actuators and possible switching 
sequences. Additionally, it would be interesting to address the sparsest feedback patterns that 
ensure the switched linear system to ensure stabilizability. 
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